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Abstract 

Temperature and luminosity functions of X-ray clusters are computed semi-analytically, combining a 
simple model for the cluster gas properties with the distribution functions of halo formation epochs 
proposed by Lacey & Cole (1993) and Kitayama & Suto (1996). In contrast to several previous approaches 
which apply the Press-Schechter mass function in a straightforward manner, our method can explicitly 
take into account the temperature and luminosity evolution of clusters. In order to make quantitative 
predictions in a specific cosmological context, we adopt cold dark matter (CDM) universes. 

Assuming the baryon density parameter f2e = 0.0125/i -2 [h is the Hubble constant in units of 
100km-sec _1 -Mpc _1 ) and the COBE normalization of matter fluctuations, temperature and luminosity 
functions of X-ray clusters depend sensitively on the density parameter f2o- Allowing for several uncer- 
tainties in observational data as well as in our simplified assumptions, we conclude that Qq ~ 0.2 — 0.5 
and h ~ 0.7 CDM models with/without the cosmological constant reproduce simultaneously the observed 
temperature and luminosity functions of X-ray clusters at redshift z ~ 0. 

Subject headings: cosmology: theory - dark matter - galaxies: clusters - galaxies: formation - X-rays: 
galaxies 
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1 Introduction 

Since clusters of galaxies are the largest virialized structure in the universe, they should retain the 
initial conditions at their formation epochs fairly faithfully. This implies that detailed studies of the 
clusters at high redshifts, as well as at present, should provide important clues to the evolution of the 
universe itself. Since X-ray identifications of clusters of galaxies are largely free from the projection 
effect which notoriously plagues optically selected cluster catalogues, X-ray observations are suitable 
for probing cosmological signatures from clusters of galaxies. Homogeneous samples of distant clusters 
of galaxies, which recent X-ray satellites such as ROSAT and ASCA are constructing, will uncover 
various statistical properties of clusters with higher reliability in the near future. Therefore, quantitative 
theoretical predictions are of great value in interpreting the observed data properly. 

Most theoretical approaches in X-ray cosmology rely on cither statc-of-thc-art numerical simulations, 
or simplified analytical formalisms. The former approach is limited by the dynamical range available on 
the present computer resources; a typical core radius of clusters is (0.1 ~ 0.2/i _1 Mpc) while the mean 
separation of the Abell clusters (richness class 1) is ~ 55ft. _1 Mpc, where h is the Hubble constant H in 
units of 100km • sec -1 • Mpc -1 . A small-scale resolution much below the core size is essential because 
a large fraction of the X-ray luminosity of clusters comes from the core. On the other hand, a large 
simulation box size is a prerequisite for statistical studies of clusters. Unfortunately, it is still hard to 
simultaneously satisfy these requirements even with the currently most advanced simulations (e.g. Kang 
ct al. 1994; Bryan et al. 1994; Cen et al. 1995). 

A major fraction of the X-ray luminosity from clusters is produced via a fairly simple process, thermal 
brcmsstrahlung. Thus one may readily compute their temperature and luminosity functions at redshift z, 
71t{T, z) and til(L, z), once the mass function um(M, z) is given, where T, L and M are the temperature, 
luminosity and mass of the clusters, respectively. Although the Press-Schechter theory (Press & Schechter 
1974, hereafter PS) is frequently applied for this purpose, it has a serious limitation in predicting the 
temperature and luminosity functions; PS theory predicts the number density of virialized objects of mass 
M collapsed before a given epoch z, but does not specify the formation epoch Zf of the objects. In fact, 
the predictions of the spherical nonlinear collapse model (e.g. Peebles 1980) suggest that the temperature 
and luminosity of objects that virialize at Z{ should scale as T(z{) oc (1 + Zf) and L(z{) oc (1 + Zf) 7 / 2 in the 
Einstein-de Sitter universe, for instance. The previous approaches based on the PS formula (e.g. Evrard 
& Henry 1991; Hanami 1993) have simply replaced Z{ by z in computing T and L (see also eq. [2.3] 
below). This procedure corresponds to implicitly assuming that T and L of individual clusters decline 
with time as T(z)/T(z f ) = (l + z)/(l + z f ) and L(z)/L(z f ) = (l+z) 7 / 2 /(l + z f ) 7 / 2 . Numerical simulations, 
on the contrary, suggest modest evolution in the opposite direction (e.g. Evrard 1990; Suginohara 1994; 
Navarro, Frenk, & White 1995). While the above assumptions may not alter n,T(T, z) and n-^(L,z) for 
larger clusters most of which would have formed only recently (zf ~ z <C 1), it is likely to affect the 
results for less massive clusters. 

This line of consideration motivates us to specify explicitly the formation epoch of virialized structures 
and their subsequent evolution in making statistical predictions for comparison with observations. A key 
quantity for this purpose is a distribution function of halo formation epochs proposed by Lacey & Cole 
(1993, hereafter LC) and Kitayama & Suto (1996, hereafter KS) in a similar but slightly different manner 
(see also Blain & Longair 1993; Sasaki 1994). In this paper, we combine these distribution functions with 
a simple model of cluster gas properties to make quantitative predictions on the temperature and lumi- 
nosity functions of clusters of galaxies in cold dark matter (CDM) universes with/without a cosmological 
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constant A . 

The plan of this paper is as follows. Section 2 outlines the formalism we use in computing the 
temperature and luminosity functions. Section 3 describes a model of X-ray clusters, and our main 
results are presented in Section 4. Finally, Section 5 summarizes our conclusions. 



2 Formulation 



2.1. Temperature and luminosity functions 



The temperature function tit(T, z) of X-ray clusters is defined as the differential comoving number 
density of clusters of temperature T at a given redshift z (the luminosity and mass functions, ni,(L,z) 
and um{M, z), are defined in a similar manner). On the basis of the spherical collapse model (see 
Section 3.1 and Appendix A), we assume that the temperature T and the luminosity L of X-ray clusters 
are determined by the mass M, the redshift of formation Zf and the redshift of observation z; i.e. 
T = T(M, Zf, z) and L — L(M, Zf, z). Then a proper theoretical prediction for tit{T, z) and nj_,(L, z) 
requires a quantity F(M, zf, z), the number density of objects of mass M which formed at Zf and are 
observed at z, rather than simply a mass function um(M, z). This is because T and L depend not only on 
M and z but also on Zf. Once F(M, zf, z) is given, the temperature and luminosity functions are written 
respectively as 

dM 
~dT 



and 



/>oo 

n T {T,z)= / dz{F(M,z f ;z) 

J z 

f°° dM 
n h (L,z)= J dz{F(M,z f ;z) — 



(2.1) 



M=M(T,z t ,z) 



M=M(L,zt,z) 



(2.2) 



By contrast, conventional approaches simply translate the PS mass function as (e.g. Evrard & Henry 
1991; Hanami 1993) 



irr, N , . dM 

n T {T,z) = n PS (M,z) — 



M=M(T,z) 



it \ , . dM 

n L (L, z) = n PS (M, z) — 



(2.3) 



M=M(L,z) 



which correspond to assuming that each cluster forms when it is observed (zi = z). In the above, the PS 
mass function nps(M, z) is given by 



n PS (M, z) 



2 Po S c (z) 
nMa 2 (M) 



do-(M) 



dM 



exp 



2a 2 (M) 



(2.4) 



where p (~ 2.78 x 10 11 fi /i 2 M Q Mpc -3 ) is the mean comoving density of the universe, a 2 (M) is the mass 
variance of linear density fluctuations at the present epoch, and S c (z) is the critical linear overdensity 
evaluated at present for a spherical perturbation to collapse at z. In what follows, a 2 (M) and S c (z) are 
computed according to the formulae presented in Appendices A and B. 

Since the conventional PS approach (eq. [2.3]) identifies the epoch of formation Zf with that of ob- 
servation z, the temperature and luminosity evolution of individual clusters after their formation are not 
properly taken into account. Our method (eqs [2.1] and [2.2]), on the other hand, can in principle include 
the evolution more naturally. For this purpose, one needs an appropriate expression for F(M, zf, z), as 
well as an evolution model which specifies T — T(M,Zf,z) and L = L(M, Zf,z). We will discuss these 
points in §2.2 and §3.1 below. 
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2.2. Distribution function of halo formation epochs 

In the hierarchical clustering scenario, each halo increases its mass via major mergers and steady accre- 
tion, and thus the formation epoch Zf of a halo is not always well-defined. LC proposed a differential 
distribution function of halo formation epochs, dp/dzf, defined via the probability that a halo of mass 
M which exists at z has a mass greater than M/2 for the first time at Zf (see LC §2.5.2): 

= (2.5) 

ozf au>{ ozf 

dp, . 1 f 1 ~ M 1 / Cj 2 \ ( u 2 \ , . 

^r(M,u> { ) =^= / dS T777 ^ (1-4 )exp(-^4) , (2.6) 

where 



dui y ' V2^Jo M'(S)S 3 / 2 \ SJ V 2S* 



UJ{(M,Zf,z) = f7=Fr f== , (2-7) 

y/a 2 (M/2) - a (M) 

5(M ' M) = a 2 (M/2)-a 2 (MY ( ) 

and M'(S) is computed by solving equation (2.8) for M'. Figure |^ shows dp/dz{ as a function of Zf in 
CDM models with Ao = 0, h = 0.7 and 6=1, where 6 is the bias parameter defined by b = l/tr(8ft.~ 1 Mpc). 

Note that, in the above definition of the halo formation epochs, the mass of a halo at Zf and z is 
different at most by a factor of 2. Supposedly this factor will be close to 2, since the increase of halo 
masses is expected to be dominated by steady accretion of small objects (although the mass increase 
is not necessarily continuous because of major mergers, the fraction of such events should be relatively 
small). Thus the LC proposal implies that the quantity F(M, Zf,z) in equations (2.1) and (2.2) can be 
written as 

F LC (M, zf, z)dMdzi = 2— -<2M, z u z)n PS (2M, z)dMdz u (2.9) 

OZf 

where M is the halo mass at Zf, which is assumed to have doubled by z, and npg(M, z) is the PS mass 
function given by equation (2.4). 

As discussed by LC, however, the probability dp/dujf given by equation (2.6) becomes negative for 
small Q{ in the case of power-law matter fluctuation spectra P(k) oc k n with index n > 0, implying 
that the above definition of the halo formation distribution function is not completely self-consistent. 
Since the effective power index of the observed fluctuation spectrum is negative below the cluster scale 
(e.g. Peacock & Dodds 1994), this would not be a serious problem for most astrophysically interesting 
objects. Nevertheless this clearly exhibits the importance of exploring other prescriptions for computing 
the distribution of halo formation epochs. 

One such alternative is given by KS. They derived expressions for the rates of formation and destruc- 
tion of bound virializcd objects based upon the conditional probability argument developed by Bond et al. 
(1991), Bower (1991), and LC. Their formalism yields the number density of haloes of mass M ~ M + dM 
that form via major mergers at Zf ~ Zf + dzf and remain in the range M ~ 2M at a later time z (see KS 
§2.4): 



Fks(M, zf, z)dMdz{ 



S c (z { ) - 25 c (z) 



x exp 



2tt yJa 2 (M/2) -a 2 {M) \ \ S c (z) 



a 2 (M) 
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rfc[Y(M,z u z)} 



d5 c (z { ) 



dzf 

where erfc(w) is the complimentary error function defined by 



n PS (M,zi)dMdzf, (2.10) 



2 f°° 

erfc(w) = -= / cxp(-i 2 )di, (2.11) 
and X(M, Z{, z) and F(Af, Zf , z) are respectively 

x(m, *, ,) , + (2 . 12) 

i/2<7 2 (M)<7 2 (2M)[<7 2 (M) - a 2 {2M)\ 

Y(M Zi z) = ajh^M (2 13) 

v /2a 2 (M)cr 2 (2A/)[a 2 (Af) - a 2 (2Af)] ' 

As mentioned by KS, the above definition (eq. [2.10]) for the formation epoch distribution may lead to a 

systematic overestimation of the number of haloes that form by mergers at Zf (see §4 of KS). 

Strictly speaking, therefore, neither _Flc nor -Pks is a completely satisfactory expression for F(M, zf, z). 

In practice, however, F^c and -Fks are in good agreement except at Zf ~ z (Fig. ^); the small discrepancy 

is ascribed to different criteria of "formation" in the two approaches, because the formation epochs are 

not uniquely defined for haloes which increase their mass quiescently via accretion. Furthermore, as will 

be shown in §4, such a difference between the LC and KS models does not affect the predictions for 

nr(T, z) and n^(L, z) (see Figs ^| and ^). Thus we believe that they provide reasonable approximations 

to F(M, Zf ; z), at least for our present purpose. 

3 Descriptions of X-ray clusters 

3.1. A simple model 

Given the distribution function of halo formation epochs F(M, Zf \ z), it remains to specify T — T{M, Zf, z) 
and L — L(M, Zf, z). We basically follow our previous model (§3.1.2 of KS), but slightly modify the 
hypothesis on the temperature and luminosity evolution. 

We make the approximation that X-ray clusters consist of dissipationless dark matter and baryonic 
gas. The dark matter component dominates the total gravitating mass and is supposed to be in virial 
equilibrium. Once the cosmological parameters such as S7q and Ao are specified, the virial radius r v ; r and 
the virial temperature T v i r are determined by the mass M and the formation epoch Zf using the spherical 
collapse model (e.g. Peebles 1980; Appendix A) We compute these quantities from the mean density p vlY 
given in Appendix A by 

rw(M,z f ) = [J- , (3.1) 

T v ,(M,z f ) = ^fmp, (3.2) 

where G is the gravitational constant, fee is the Boltzmann constant, m p is the proton mass, and [i is the 
mean molecular weight. Hereafter we assume that the intracluster gas is fully ionized with primordial 
abundances of helium and hydrogen, and thus set /i = 0.59. 

For a cluster temperature of > 3keV, the contribution from the line emission on the total X-ray 
luminosity can be neglected. Since the shape and amplitude of the temperature and luminosity functions 
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are derived observationally for relatively high temperature clusters (see the error box in Figs ^| ~ [T^ below; 
Henry & Arnaud 1991), we compute the X-ray luminosity only from thermal bremsstrahlung emission. 
Incidentally, the X-ray line emission from less massive clusters may be important in considering the origin 
of the soft X-ray background (Cen et al. 1995; Suto et al. 1996), and we will discuss the effect of line 
emissions elsewhere (Sasaki et al. 1996). 

We assume that the cluster gas is isothermal with temperature given by 

T(M, z f , z) = k{zu z)T vir (M, z { ), (3.3) 

where n(z{ , z) is introduced so that temperature evolution due to quiescent accretion of matter after Zf can 
be taken into account. Note that this treatment is possible because our method explicitly distinguishes 
Zf from z. As there is not yet a good theoretical model to describe the temperature evolution, we simply 
suppose that n(zf , z) takes a power-law form: 

K ( zuz) =(i±^y . (3.4) 

Hydrodynamical simulations indicate that the temperature of clusters, once formed, does not change 
drastically after virialization, and are roughly consistent with < s < 1 (e.g. Evrard 1990; Navarro, 
Frenk & White 1995). 

We adopt the following spherically symmetric distribution for intracluster gas density (e.g. Jones & 
Forman 1984) 



PgasM = Pg 



ga.s 



1+1 - 



2 

r ' 



(3.5) 



where /9g as is the central gas density, and r c is the core radius. Since this is basically an empirical fitting 
formula, it is difficult to predict r c . Thus we adopt a simple self-similar model in which r c is proportional 
to r vir (M, z f ): 



r c = O.lS^Mpc x J:^ lZf) n , (3.6) 

r vir (lO 15 M ,z f = 0) V 1 

where the normalization is chosen to match the X-ray observations (Abramopoulos & Ku 1983; Jones & 
Forman 1984). Now /9g as is fixed by the relation: 

P gas (r)^r 2 dr = M(^j , (3.7) 

where J7b is the baryonic density parameter of the universe chosen to be consistent with primordial 
nucleosynthesis: £!b = 0.0125/i~ 2 (Walker et al. 1991). We assume no intrinsic evolution for the gas 
density profile, which is also suggested by numerical simulations (e.g. Navarro, Frenk & White 1995). 

Combining all the above assumptions, the temperature T and the bolometric luminosity Lboi of X-ray 
clusters in our model are roughly related to M , Zf and z as follows: 

T(M, z f , z) oc M 2 / 3 (l + z f )« (3.8) 

L bol (Af,Zf,z)cxAf 4 / 3 (l + z f ) 7 «/ 2 (i±^y /2 , (3.9) 

where £ is the effective index which varies weakly with Zf, Qq and Ao; e.g. £ = 1 independent of Zf in the 
case of (fi 0) A ) = (1,0), while f ~ 0.5 at z f - and f ~ 1 at z f > 2 for (f2 ,A ) = (0.2,0.8). 



n 
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3.2. Comparison with observed correlations 

The model presented in the previous subsection is admittedly crude, so we need to examine the extent 
to which it is consistent with the available observations before presenting the results for n^(T^z) and 
nt,(L, z). Specifically, we plot the correlations among T, r c and L(2 — lOkeV) in Figures ||, || and [| In 
these figures we select three representative CDM models, Eh5, L3h7, and 03h7, which have respectively 
(f2o) Ao, h) = (1.0, 0, 0.5), (0.3, 0.7, 0.7), and (0.3, 0, 0.7). In fact, model L3h7 is shown to be reasonably 
consistent with the observed temperature and luminosity functions at z — in the next section. 

Basically, our models reproduce fairly well the observed T — r c and L{2 — lOkeV) — r c correlations, 
considering that most clusters present at z = would have formed at Zf < 1. The slopes of our models in 
Figures |] and |I] would become less steep if the distribution of Zf is taken into account, because massive 
clusters (large T and L) are likely to have formed more recently than less massive ones. As regards 
the L(2 — lOkeV) — T correlation on the other hand, only low f2o models show reasonable agreement 
with the observed data at large L and T . The predicted curves are less steep than the observed one. 
In fact, several authors have pointed out that simple adiabatic models fail to account for the observed 
L — T relation (e.g. Navarro, Frenk & White 1995). It has also been suggested that the discrepancy may 
possibly be reconciled if the gas has attained large entropy at high redshift (Kaiser 1991; Evrard 1990; 
Evrard & Henry 1991). 

Apparently this is an important area of research, especially if one takes seriously the above discrep- 
ancy between the observed L — T relation and the theoretical prediction of a simple adiabatic model. 
Nevertheless it is fairly independent of our current formalism; once any physical model of such pre- heating 
specifies T — T(M, Zf, z) and L = L{M, Zf, z), then we can compute the resulting tit(T, z) and nh{L, z) 
by a straightforward application of the method described here. Thus we do not discuss this problem 
further and hope to come back to it elsewhere. 



4 Predictions for temperature and luminosity functions 

The temperature and luminosity functions of clusters predicted in CDM universes at z — are shown in 
Figures ^ and 0, respectively. All the models adopt h = 0.7 based on the recent observations of Cepheids 
by the Hubble Space Telescope (Tanvir et al. 1995). The fluctuation amplitudes are normalized according 
to the COBE 2 year data (Gorski et al. 1995; Sugiyama 1995). Panels (a) and (b) in these figures are based 
on the LC model (eq. [2.9]) while panels (c) and (d) adopt the KS model (eq. [2.10]). The comparison 
implies that the results are quite insensitive to the choice of models, because the small difference between 
the two models shown in Figure [2] has been canceled out in the integration over zt (eqs [2.1] [2.2]). Thus in 
what follows we will adopt the LC model just for definiteness. Clearly the predictions of the temperature 
and luminosity functions vary largely with the shape and the amplitude of the fluctuation spectrum. 
Since both of these are fixed by £lo and Ao assuming the CDM model and the COBE normalization, 
comparison of the theoretical predictions with the observed temperature and luminosity functions serves 
as a possible discriminator of the values of FIq and Ao- 

We have tested the robustness of the above predictions against our assumptions about the temperature 
evolution in Figure^. The temperature functions depend rather sensitively on the evolution index s. In 
fact, the temperature evolution assumed in our model, T(z)/T(z{) = (1 + Zf) s /(1 + z) s with s = 1 (thick 
lines), tends to move the predictions of ut(T 7 z = 0) to the right roughly by a factor of (1 + zi) relative 
to the s — case (thin lines). Although the formation epoch Z{ varies with M and therefore with T of 
an individual cluster, (zf) ~ 0.5 for typical clusters reasonably accounts for the behavior in Figure ||a. 
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On the other hand, the shift in the luminosity function is rather small (Figure |]b), because our evolution 
model assumes that the cluster density is fairly constant after virialization and the bolometric luminosity 
evolves as L hol {z) / L hol (z f ) = (1 + z f ) s / 2 /(l + z) s / 2 . 

Incidentally, Figure || shows that our s = model prediction is very close to the original PS approach 
(eq. [2.3]), especially for the temperature function. However, this does not imply that the effect of zi ^ z 
is negligible. Our predictions have been shifted relative to the original PS approach in the following two 
ways. One is the horizontal shift due to the different hypotheses for the temperature and luminosity 
evolution; the original PS approach implicitly assumes for instance T(z)/T(z{) = (1 + z)/(l + Zf) and 
L(z) I L(zf) = (1 + z) 7 / 2 /(l + Zf) 1 / 2 in the f2o = 1 universe, while we assume those just mentioned in 
the last paragraph. The other is the vertical shift due to the different mass of haloes at Zf and at z; our 
formalism assumes that the mass of a halo at Z{ is half of that at z, while no mass change is considered 
in the original PS approach. The former shift depends on the evolution index s, and the latter one is 
sensitive to the fluctuation spectral index. On cluster scales, the above effects have almost canceled each 
other for the temperature function with s = 0. 

In any case, Figures |[ fj] and || indicate that, with the COBE normalization, a change in the parameter 
s in the range < s < 1 roughly corresponds to changing Slo by 0.1 ~ 0.2. In the case of s — and 
h = 0.7, either the (f2o, Ao) — (0.4, 0.6) or the (0.5, 0) CDM model fit best within the observational error 
box. If s = 1 is adopted instead, models with (f2o, Ao) — (0.3,0.7) or (0.4,0) are prefered. Of course, the 
conventional argument of the age problem also points to a lower Qq; if h — 0.7 as we adopted here, the 
ages of the universe are 13.5 Gyr, 12.4 Gyr, 10.9 Gyr and 10.5 Gyr for (O ,A ) = (0.3,0.7), (0.4,0.6), 
(0.4, 0) and (0.5, 0) respectively. 

Adopting the COBE normalization in the computation of the temperature and luminosity functions, 
we implicitly assume that luminous clusters trace the underlying mass distribution. Although this seems 
quite reasonable on scales relevant for formation and clustering of galaxy clusters, it is not easy to 
rigorously justify this assumption. So we change the normalization of the fluctuations and plot the 
results for b = 1.5 and 0.5 in Figure ^ (assuming s = 0.5). For comparison, the COBE 2 year data imply 
that b = 0.73, 0.87, and 1.8 in models Eh5, L3h7 and 03h7 respectively. 

With the above uncertainty and the observational error in mind, an acceptable range of Oo from 
the observed temperature and luminosity functions is 0.2 Slo ~ 0.4 in the case of spatially flat CDM 
universe, and 0.3 flo ~ 0.5 for open CDM universe. It should be noted that the Slo = 1 CDM model fails 
to account for the observation by a wide margin and can be ruled out even from the present comparison 
only. 

More interesting for the future X-ray observations are the predictions of evolutionary behaviors in the 
temperature and luminosity functions. Figure [II] plots the predictions at z = and 1. Evolution of the 
temperature function (Figure |l0|a) is in general negative (i.e. decreasing number density toward the past) 
for massive clusters (large T), while it becomes positive for smaller ones. In model L3h7, for instance, 
a significant negative evolution occurs at T > 2keV, where observed data are available at z ~ 0. On 
the other hand, the luminosity function mostly shows very weak evolution in the range of observational 
interest (Figure [To|b) . This is still consistent, within the observational error bar, with the results by 
Gioia et al. (1990) which indicate weakly negative evolution at L(0.3 — 3.5keV) > 10 44 /i _2 erg s _1 . Our 
model can be tested against future observations through the different evolutionary signature of h.t(T, z) 
(negative evolution) and 71l(£, z) (almost no evolution). 
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5 Conclusions 

We have examined statistical properties of X-ray clusters in CDM universes in a semi-analytic manner, 
which is complementary to numerical simulations. Our method is different from the previous approaches 
based upon the conventional PS theory (e.g. Evrard & Henry; Hanami 1993), in that we explicitly took 
account of the epochs of cluster formation adopting the halo formation epoch distributions proposed by 
LC and KS. As a result, our method can include the subsequent evolution of cluster temperature and 
luminosity. 

Although the LC and KS proposals involve slightly different definitions of the halo formation epochs, 
we found that the resultant temperature and luminosity functions are remarkably similar. Deviations 
from the previous PS approach become larger for clusters that form at higher redshift and for stronger 
evolution. 

The shape and amplitude of the temperature and luminosity functions vary sensitively with the cos- 
mological density parameter Q , if the fluctuation amplitude is fixed by the COBE normalization. Given 
the qualitative nature of our simple model, however, we should not constrain the cosmological parameters 
so stringently. Rather we conclude that the low-density fio ~ 0.2 — 0.5 CDM models with/without the 
cosmological constant are roughly consistent with the observed temperature and luminosity functions. 
Nevertheless we can argue against Oq = 1 CDM models at least from the present analysis only. 



We thank Takahiro T. Nakamura and Naoshi Sugiyama for useful discussions on the spherical collapse 
model and the COBE normalization. We arc grateful to Ewan D. Stewart for a careful reading of the 
manuscript. This research was supported in part by the Grants-in-Aid by the Ministry of Education, 
Science and Culture of Japan (07740183, 07CE2002). 
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APPENDICES 
A Spherical collapse in fi < 1 universes 

In the present analysis, we need the mean density p v i T (z{) and the critical linear overdensity S c (zf) of 
an object that virializes at Zf in a universe with arbitrary fio (fi < 1). For definitencss, we summarize 
the results of the spherical collapse model as well as the fitting formula we used for Ao = 1 — fin models 
(Peebles 1980; Lahav et al. 1991; LC; Nakamura 1996). Note that S c given below is the value computed 
at virialization, and 5 c (zf) used in our main text is the value computed at present. They are related via 
S c (z{) = S c ■ D(z{ = 0)/D(z{), where D(z) is the linear growth rate. 



1. fin = 1, 



2. fi < 1, A = , 



Pvir(zf) 
Sc 



= 18tt 2 ~ 178, 

3(12^) 2/3 
20 



Pvii(zf) _ ^ 2 (cosh7?f - l) 3 
p(z f ) (sinh77 f - 77 f ) 2 ' 

3 sinh?7f (sinh?7f — rjf) 



(coshr^f — l) 2 



1.69, 



2tt 



sinh rjf — r]{ 



2/3 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



3. fin < 1, An — 1 — fin > 



PvirCgf) 



rta \ 2w f 
18^ 2 (l + 0.4093 Wf a9052 ), 



Sc = 



l F (Wi'- W{ 



2^y/ 3 

X J 



3fl27r1 2 / 3 
1 2Q > (1 + 0.123 log 10 fi f ), 



(A.5) 



(A.6) 



In the above, p(z{) = po(l + Z{) , T}{ = arccosh(2/fif — 1), Wf = 1/fif — 1, x — ^oHo r ta/ \GM) , r ta is the 
maximum turn-around radius and F is the hyper-geometric function of type (2,1). The density parameter 
fif is defined at Zf and is related to fin and Aq through 



fi f = 



fi (l + z f ) 3 



fi (l + z f f + (1 - fi - A )(l + zf) 2 + A ' 



(A.7) 



B Mass variance for the CDM fluctuation spectrum 

The mass variance cr 2 (M) is related to the linear power spectrum of density fluctuations P(k) through 

1 r°° 

<j 2 {M) = / P(k)W 2 (kr M )4nk 2 dk, (B.l) 

( 27T r Jo 
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where the filtering radius r M is related to the mass M by r M — [3M/(47rp )] 1 / 3 as the top-hap window 
function is used: 

3 

W(x) = — r(sin:z; — xcosx). (B.2) 

x 6 

We adopt the CDM power spectrum of the form given by Bardeen et al. (1986) with the scale invariant 
initial condition: 

2 

[1 + 3.89<7 + (16.1g) 2 + (5.46<?) 3 + (6.71g) 4 ] ~ 1/2 , (B.3) 

where q = k/{Th Mpc -1 ). Here the quantity V for non-negligible baryon density is given by (Peacock & 
Dodds 1994; Sugiyama 1995) 

T = Q h(T a /2.7 K)- 2 exp[-0 B (l + V2MIq % (B.4) 

where To is the temperature of the cosmic microwave background radiation. 

We found that a 2 and da 2 /dM calculated from the above equations are accurately fitted simultane- 
ously by the following formula and its derivative with respect to m: 

a 2 oc [1 + 2.208m p - 0.7668m 2p + 0.7949m 3p ] , (B.5) 

where p = 0.0873, and m = M(Th) 2 / (10 12 Mq). The above approximation holds within a few percent 
for both a 2 and da 2 /dM in the range 10~ 6 < m < 10 4 . Throughout the present paper, a is evaluated 
via equation (B.5) and normalized by 

er(r M = M^Mpc) - 7, (B.6) 


where b is the bias parameter. 



P(fe) oc k 



ln(l + 2.34g) 
2.Mq 



C Fitting formulae for the Lacey & Cole distribution function 



Since the distribution function of the halo formation epochs derived by LC (eq. [2.6]) requires a time- 
consuming numerical integration, the following fitting formula is used in the present paper. 

First we consider the case of power-law fluctuation spectra P(k) oc k n , i.e. a 2 (M) oc M~ a where 
a = (n + 3)/3. Then equation (2.6) reduces to 



9p, ~. /8-o?/2 2(2" -1) , 



/±± y 

V2 



1 



(C.l) 



For u>{ <C 1, equation (C.l) becomes 

dp , . . 



V 7T a Jo 



2(2 a -l)„ 
+ - -a) f + OK 2 ), 



(C.2) 



For a)f 3> 1, on the other hand, we found that numerical results for equation (C.l) are insensitive to a 
and well approximated by the solution in the a = 1 case: 



TT-(a — l,<Df) = 2o)f erfc 



(C.3) 
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On the basis of equations (C.2) and (C.3), we constructed the following interpolation formula which 



gives a very good fit to equation (C.l): 

dp A(a) 

acjf 1 + -B(a)cL>f 

where 



2C(a)£ f erfc ( -j= 



A{a) 
B{a) 
C(a) = 1 



-(1 -a)(0.0107 + 0.0163a), 

7T 



A(a) 



C(a) 



2" - 1 



1-a 
25 



(C.4) 

(C.5) 
(C.6) 
(C.7) 



Since equations (C.1)-(C.3) are valid only when the underlying spectra are of power-law form, the 
above fitting formulae may not be applicable to more general spectra. In fact, in the case of the CDM 
spectrum, a straightforward substitution of the effective spectral index a c s = — c?lncrQ DM /dlnM into 
equation (C.4) does not work very well, because a e g varies with the mass scale and so cannot be approx- 
imated by a constant in the integral of equation (2.6). Nevertheless, we found that equation (C.4) still 
gives an excellent empirical fit to equation (2.6) for the CDM spectrum, if <S ff below is adopted rather 
than a e g: 

& cS = a GfT (0.6268 + 0.3058a ofr ). (C.8) 

The accuracy of this approximation is within a few percent in the range 10 6 M Q < M(iloh 2 ) 2 < 10 16 M 
(see Fig. P). 
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Figure 1: The formation epoch distribution function dp/dzf by LC in CDM models; (a) Clo = 1, z = 2, 
(b) fi = 1, z = 0, (c) fi = 0.1, = 2, (d) fl = 0.1, z = 0. In all panels, A = 0, /i = 0.7, and b = 1. 
Lines indicate M = 10 9 M o (solid), 10 12 M Q (dotted), and 1O 15 M (dashed). 



KlTAYAMA & SUTO 



15 



N 

< 



O 
Q. 



O 



3 



O 



i i i — i — i — r 



AAAAAAAAAAA^ 




Figure 2: The distribution function of halo formation epochs by LC (eq. [2.9]) and KS (eq. [2.10]) in 
CDM models; (a) O = 1, z = 2, (b) Q = 1, z = 0, (c) ft = 0.1, = 2, (d) Q a = 0.1, = 0. In all 
panels, A = 0, h = 0.7, 6=1. Lines and symbols indicate the LC and KS formulae respectively for 
M = 10 9 M Q (solid line, triangles), 10 12 M Q (dotted line, crosses), and 10 15 M (dashed line, circles). 



1G 



STATISTICAL PROPERTIES OF X-RAY CLUSTERS 




r c [h 1 Mpc] r [h 1 Mpc] 

Figure 3: The T — r c relation for the clusters in CDM models at (a) z = and (b) z — 0.5. Three 
representative models, Eh5, L3h7, and 03h7 are plotted in solid, dotted, and dashed lines, respectively. 
Thick and thin lines correspond to the results for Z{ = z + 1 and Zf = z, respectively. The crosses show the 
Einstein observations (Jones & Forman 1984; David et al. 1993), and the diamonds indicate EXOSAT 
data (Edge & Stewart 1991a,b). 
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Figure 4: Same as Figure except for the L(2 — lOkeV) — r c relation. 
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Figure 5: Same as Figure || except for the L(2 — lOkeV) — T relation. 
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T [keV] T [keV] 

Figure 6: Predictions of the temperature function of clusters at z = in CDM universes with h = 0.7; 
(a) LC model with A = 1 — Sl , (b) LC model with A = 0, (c) KS model with A = 1 — flo, and (d) 
KS model with A = 0. The fluctuation amplitudes are normalized according to the COBE 2 year data 
(Gorski et al. 1995; Sugiyama 1995). The thick solid box indicates the observational error box (Henry & 
Arnaud 1991). 
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Figure 7: Same as Figure except for the bolometric luminosity function. 
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Figure 8: Effect of different evolutionary assumptions on the predictions of (a) temperature and (b) 
luminosity functions at z = in Eh5, L3h7, and 03h7 CDM models (COBE normalized). Results for the 
s = 1 and s = models are respectively plotted in thick and thin lines, while those based on the original 
PS model are plotted in symbols. 




Figure 9: Effect of fluctuation amplitudes on predictions of (a) temperature and (b) luminosity functions 
at z = in Eh5, L3h7, and 03h7 CDM models (s = 0.5). Thick and thin lines indicate the results for 
b = 1.5 and b = 0.5, respectively. 
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Figure 10: Evolution of (a) temperature and (b) luminosity functions in Eh5, L3h7, and 03h7 CDM 
models (s = 0.5). Thick and thin lines indicate the results for z = 1 and z = 0, respectively. 




Figure 11: The scaled distribution function of halo formation epochs, dp/dujf, in the standard CDM 
model (fi = 1, Ao = and h = 0.5); (a) linear and (b) log-log plots. Symbols and lines indicate results 
of a numerical integration of equation (2.6) and a fit by equations (C.4)-(C.8), respectively; M — 10 9 Mq 
(triangles, solid line), 10 12 M (crosses, dotted line), and 10 15 M Q (circles, dashed line). 



